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Abstract: We compute a shift of baryon mass spectra due to quark masses in 
perturbation, in Sakai-Sugimoto model of holographic QCD. We find the shift for 
the ground state nucleons to be SM = 4.1 GeV~^, which is consistent with the 
current lattice QCD result. We predict the same value of the shift for N(1535) and 
A, while a larger value 7.7 m^GeV"^ for Roper N(1440). We also present some 
evidences that the shifts of the vector meson masses are suppressed in the large 't 
Hooft coupling limit. 
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1. Introduction 

Recent progress in holographic QCD of string theory enables us to compute the im- 
portant physical observables of QCD: the hadron masses. QCD has only two kinds of 
dimension-ful parameters: the prime one is the QCD scale which basically generates 
all the hadron masses, and the quark masses. For the light quarks, the quark masses 
can be treated as a perturbation (while heavy quarks have their particular impor- 
tance in various phases in QCD). In this short paper, we present a computation of 
the hadron mass shift due to the quark masses, in the Sakai-Sugimoto model 0, 
which is the most successful model of holographic QCD so far. 

Sakai-Sugimoto model is a holographic dual of massless QCD. There proposed 
two ways to introduce the quark mass to the model: (i) worldsheet instantons Q 
(see also [§]), and (ii) tachyon condensation of the D8/D8-branes 0. We are going 
to use the worldsheet instanton method (i), since the other method (ii) assumes a 
tachyon "effective" action which is difficult to be validated in any manner in string 
theory except for using string field theories.* 

The worldsheet instanton produces the quark mass term in QCD, and the gravity 
dual of the worldsheet instanton generates a term of the form 

6L = ctT[MU] (1.1) 

*The tachyon approach has more direct analogue with a famous bottom-up model of holographic 
QCD and the mechanism is similar to the Higgs mechanism, while the worldsheet instanton 
approach mimics extended technicolor, see [|[ ^. 
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where M is the quark mass matrix of the form M =diag(mu, Tn^, and U is the pion 
field in the standard notation in the chiral perturbation theory: U = exp [2i7r (x) //^] . 
The value c is computed from the worldsheet instanton amplitude, but it turned out 
to be difficult to evaluate (for some computations, see @, H, ^). Except for the pion 
mass through (|1 . 1| ) , the contribution of quark masses to other meson masses has not 
been computed yet. 

In this paper, we present two computations. First, in Sec. |^, we compute shift 
of the baryon mass spectra due to the quark mass, by assuming that the quark mass 
dependence in the meson effective lagrangian of the Sakai-Sugimoto model appears 
only in the new term (|1 . 1|) , at the leading order in expansion in 1/A. This assumption 
will be discussed in Sec. ||. There we will provide some computations and show that 
shift of the vector /axial vector meson masses is small and at a higher order in 1/A 
expansion. 

Under the assumption mentioned above, it is quite easy to compute the baryon 
mass shift. As noted in and computed in detail in [|10[;^ the baryon in the Sakai- 
Sugimoto model is nothing but an instanton-like soliton in the 5-dimensional Yang- 
Mills- Chern- Simons (YMCS) action describing all the mesons effectively in a unified 
way. The term ( p..lD can be thought of as a perturbation to the YMCS action, so 
the mass shift can be computed by just inserting the soliton configuration obtained 
in [jlO| to ( |1 . 1|) . The resultant correction is evaluated with baryon wave functions of 
[0. Since the parameter c is undetermined, we obtain the baryon mass as a function 
of the pion mass. Our result for nucleons is in rough agreement with results of lattice 
QCD. The mass shifts for N(1535) and A turn out to be the same as that of the 
nucleons, while the shift for Roper N(1440) is larger. 



2. Baryon Mass Shift 

2.1 Set-up in Sakai-Sugimoto model 
2.1.1 A review of baryons 

The Sakai-Sugimoto model |l|, |^ is described by the following five- dimensional U (Nf) 
YMCS action in a curved background, 

" 1 



S = Sym + Scs , Sym ^ J ^-'^^^^ tr 



2 



[2.1) 



Scs is the 5-form Chern-Simons term 5'cs- The U{Nf) gauge fields A^{fi = 0, 1, 2, 3) 
and Az have their field strength defined in the standard manner, JF^^, = df^iAu — 
duAfj_—i[Afj_, Au]- The extra dimension z unifies all the vector /axial vector mesons via 
Kaluza-Klein like decomposition of the gauge fields, through the nontrivial functions 



^See also [Til , |l2| for an alternative description of baryons by introducing 5-diinensional fermion 
fields. Other work for introducing baryons in holographic QCD includes |13[. 
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h{z) = (1 + 2;^)^^/^ and k{z) = 1 + z^. The constant k is defined as k = aXN^ 
with a = l/(2167r'^), where A is the 't Hooft couphng constant. The dimensionful 
parameter Mkk in the action (|2.1|) is already put to be the unity, and it was argued 
in [0, that the values Mkk = 949 MeV and k = 0.00745 can fit the observed values 
of the p meson mass rrip ^ 776MeV and the pion decay constant ^ 92.6 MeV. 
In this paper we consider Nf = 2, the two-flavor case, and then the U{2) hermitian 
matrix A is decomposed by Pauli matrices t"' [a = 1,2,3) and a unit matrix I2 = r°. 



A = A + A^ 



T°- - I2 

A''— + A — 
2 2 



(2.2) 



The baryon is identified as a soliton solution localized in the 4-dimensional 
space (M = 1,2,3,2;) [^. The instanton number of the Yang-Mills theory ( p^.l] ) is 
identified with the baryon number. 

The authors of found an explicit solution of the equations of motion of ( |2.1| ). 
In view of the fact that the size of the baryon scales as A~^/^ as found in |Tl], p!0|] , 
an appropriate rescaling of the gauge fields and the coordinate by a factor of A was 
identified in |T^. But here, for our later purpose, we write the un-rescaled solution. 
The solution is almost identical to a BPST instanton except for a 1/A "correction" 
appearing in the U{1) part of the zeroth component of the gauge field. This sub- 
leading contribution can be understood from the fact that the solution represents a 
baryon whose charge should be identified with the electric charge of the trace of the 
vector part of the chiral symmetry. The explicit solution is 



cl 



1 1 



and the non-Abelian part is given by the BPST instanton. 



Ao = Am = 0, (2.3) 



/(O 



Z)l2 + i{x - X) ■ T 



[2A) 



Here the moduli parameters of the instanton come in explicitly: the instanton size 
p and the instanton location = {X'^ , X'^ , X^ , Z) = {X,Z). ^ is the distance in 

this 4-dimensional space, ^ = \J [z — Zy + \x — Xp. In addition to these 5 moduli 
parameters, the rotation in the gauge SU (2) space are moduli, which are incorporated 
as 1/ e SU(2), 



A 



M 



VA%V-^ - iVSmV 



(2.5) 



Using the standard technique for quantizing solitons |1^ , quantized baryon states 
appear in the degrees of freedom of {Z,p, V), a part of which acquires definite po- 
tentials in their quantum mechanics via the effect of the curved spacetime and S'cs 
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[0 . For example, the classical size and the location Z of the instanton can be fixed 
by minimizing the potential as 

1 1 /6 



Pci = i^\/^, ^ci = 0, (2.6) 

so the size of the baryon is (9(A~^/^). The spectrum of the baryon can be computed 
by evaluating the hamiltonian of (Z, p, V^) with their quantum wave functions. 

2.1.2 A review of quark mass 

We introduce the quark mass to the Sakai-Sugimoto model by the worldsheet in- 
stantons. The method was first introduced in and In these papers, there are 
some differences in regularization of the worldsheet instantons, and in this paper we 
follow i. 

The quark mass can be introduced by joining the D8-brane and the D8-brane by 
D6-branes Q. Once the D8 and the D8 are joined in this way, it is possible to put a 
Euclidean string worldsheet whose boundary is defined by the color D4-branes, the 
D8-brane, the D6-brane and the D8-brane: these four sets of D-branes form a square, 
in the middle of which one can put the worldsheet instanton. Since the left-handed 
and right-handed quarks live at the D8-D4 and D8-D4 intersections respectively, the 
worldsheet instanton involves with qhlR vertex which is nothing but the quark mass 
operator. The instanton amplitude, which is proportional to exp[— y4/27rQ;'] where A 
is the area of the square, gives the quark mass. Therefore the quark mass can be 
chosen by tuning the location of the D6-branes and their number. 

In the gravity dual, the introduced D6-brane is still a probe D-brane, so one can 
put the worldsheet instanton in the same manner, but now two of the corners are 
smeared out by the background curved geometry. The meson excitation described 
by the gauge fields on the D8-branes induces a worldsheet boundary coupling in the 
worldsheet instanton amplitude. 



6S = c d^x Ptr 



M exp 



A-Az 



+ C.C., M^("^" ° )(2.7) 
' met 



where and are quark masses, and z = ±Zm is the location of the D6-brane 
{z = Zm and z = —z^ specify the same radial position, on the D8 and the D8). 
This exponential form is a familiar boundary coupling of string worldsheet, and the 
integral is over the period where the string worldsheet ends on the D8-brane. For 
nonzero values of the quark mass, the chiral symmetry is explicitly broken, which is 
consistent with the standard expectation for the quark mass term in QCD. In a 
rough evaluation of the constant c was given. The constant subtraction — 12 is for 
our convenience, to make sure that when Az = the instanton amplitude vanishes. 
The Az field in nothing but the pion field H , and the relation is 



Pexp 



00 



~i / dzAz 



exp [2tn{x)/U] = U , (2.8) 
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where U is the standard notation for the pion field in the chiral perturbation theory. 
Using this, the additional action to the Sakai-Sugimoto model is 

5S = j d^x 5L , 5L = ctT [M{U + - 2I2)] , (2.9) 

which is the well-known term in the chiral perturbation theory. Here although is 
a finite value, we approximate the exponential part in ( p.7|) by ( p.8|) . For the purpose 
of a rough estimate, this simplification can be justified since the pion wavefunction 
is localized at 2; = 0. 

Using the pion mass m^, we obtain a relation 

c(m„ + m,) = ^f^ml. (2.10) 

We are working in the normalization 7i{x) = TT°'{x)Ta, tr[TaTb] = ^6ab- 

We have seen that the introduction of the D6-brane ending on the D8 and the 
D8-branes can induce a quark mass term. However, the additional action ( |2.7] ) should 
not be the only correction to the Sakai-Sugimoto action. In fact there are possibly 
two other sources for new terms coming from the presence of the D6-brane: 

• Vector meson / axial vector meson mass terms. 

• Shift of D8-brane configuration pulled by the worldsheet instanton. 

The first one is very interesting, as it may give an interesting result on shift of the 
meson mass spectra due to the quark mass. However, in the next section, we are 
going to discuss that this shift is suppressed by 1/A in the holographic QCD. The 
second one, the shift of the location (shape) of the D8-branes, is intriguing of its own 
and was studied in 0]. Since this is obviously a secondary effect in the sense that 
the string worldsheet instanton pulls the heavy D8-branes, we will not consider it in 
this paper. 

In the following in this section, we will ignore these two additional possibilities 
and consider only ( |2.7[ ) for computing the baryon mass shift. 

2.2 Baryon mass shift : classical evaluation 

When the quark mass is small, at the leading order in the quark mass, the shift of 
the baryon mass is simply given by 

SM = - [ d^x 5L[A''^] (2.11) 



■I- This tiny modification of the shape of the D8-brane can be numerically computed, and in fact 
it increases the mass of baryon. See for a relevant discussion. This will also change the meson 
masses. 



- 5 - 



where A'^^ is the classical solution (|2.3| ) representing the baryon, and 6L is given by 
( p.7|) . By the additional term ( p.7|) , the baryon configuration itself would be modified 
and it would shift the baryon mass at the higher order in the quark mass, but we 
neglect the higher order terms in this paper. 

If we substitute the classical solution to (|2.7|) , the resultant expression is a func- 
tion of the instanton moduli. At the classical level (which is the large limit), 
we can substitute the classically fixed moduli (|2.6|). To distinguish different baryon 



states, we use baryon wave functions of the moduli given in [rU| and take the ex- 
pectation value of ( |2.7| ), which gives the baryon mass shift for a given baryon state. 
Here we shall compute the classical value first, and then in the next subsection we 
will compute the quantized value for each baryon state. 

Let us evaluate this ( p.ll|) explicitly. First, note that in the large A expansion, 
the leading order term of the solution is just the same as the BPST instanton, while 
the term of the next-to-leading order is nonzero only for the zero-th component 
Aq{x, z) of the overall U{1) part of the gauge field. On the other hand, the integral 
in defining U in ( |2.8| ) which appear in (|2.11| ) is only for the z component. Thus, the 
sub-leading term does not contribute to our calculation, so what we need is just the 
BPST instanton part. 

Eventually, the integral ( ^78|) was already evaluated in [0, so let us use it. For 
simplicity, we put all the moduli to be the classical values, (|2.6| ) and V = I2, to obtain 
a classical shift of the baryon mass (we can put X° = without losing generality). 
We work in a singular gauge of the BPST instanton, 

1 , , 

(2.12) 



The Az configuration (|2.12|) is proportional to a matrix (x")r'* for any value of 2, so 
the evaluation of the path ordering in ( p.ll| ) reduces to an Abelian problem . We 
obtain 



U = exp [i/(r)x"r"] 



with 



/(r) = 71 



1\2 j_ (rr.2\2 j_ {rf.3\2 



With this expression, we obtain 

U + U^ = 2cos(/)l2 
Therefore, using the relation ( p.lO|) , we have 

6M= [ Sx flmlil-cosU)) 



x^r + [X 



(2.13) 



(2.14) 



(2.15) 



ATrf^ml.p^ dr { 1 — cos 



TT 1 



(2.16) 
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where we made a change of a variable as r 
evaluated as 



pr. The last integral is numerically 



dr 



cos 



TT 1 



1.104. 



(2.17) 



Using this, we obtain a formula for the shift of the baryon mass spectrum due to the 
quark mass, 



5M 



4:7; f I X 1.104 . 



mt 



(2.18) 



In the classical limit A^^^ 



oo, p is equal to pci which is given in (|2.6|) , so 

3/2 



5M 



1.104 X A'Kflmlpl^ = 1.104 x ^Trf^ml ' ^'^^^ 



(2.19) 



In this limit, the mass shift ( |2.18| ) is independent of the baryon state, since the 
dependence on the baryon state appears at the quantum level which is sub-leading 
in the 1/Nc expansion. In the next subsection we compute the l/N^ corrections and 
the baryon state dependence. It was found that 1/Nc corrections are rather large for 



various static properties of baryons (such as charge radii) in |jT8 



Note that our baryon mass shift SM is linear in as seen in ( 2.18| ), to the 
lowest order in m^. This is in agreement with the lattice QCD results and also with 
chiral perturbation theories. 

2.3 Baryon mass shift : baryon state dependence 

The classical evaluation of the baryon mass shift, ( |2.18| ), will be corrected once we 
take into account the quantum states of the baryon. Each baryon state is specified 
by the quantum numbers {/(= J), I3, Up, nz} where I is the isospin, J is the spin, Up 
and nz are quantum numbers associated with the moduli parameters p and Z [p^] . 
In general the mass shift depends on these quantum numbers. 

The dependence can be evaluated by the standard perturbation in quantum 
mechanics. The additional hamiltonian is ( p.ll|) evaluated with all the moduli de- 
pendence in Aci- 

First, let us consider Z-dependence. Since the Z-dependence comes in as z — Z 



in the solution ( |2.3| ) with (p.4|), once we perform the integration (|2.8|), Z disappears. 
So there is no additional hamiltonian for Z. 

Next, let us consider V. This SU(2) rotation V appears in the solution as a gauge 
transformation of the solution, (|2.5| ), the transformation of U is given as VUV~^, 
soU + W ^ V{U + W)V-'^. Now, looking at (|2l5| ), we notice that in fact this V 
dependence disappear, because U + W is proportional to the unit matrix. So, there 
is no additional hamiltonian for V. 
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Finally, let us consider the moduli p. The expression ( p.l8|) has a factor p'^, so 



it is the additional hamiltonian for the quantum mechanics for p. In sum, the shift 
of the baryon mass is dependent only on the quantum number Up. It is not sensitive 
to spin/isospin and also the quantum number nz- This means that, in particular, 
proton, neutron, A^(1535) and delta excitation have the same mass shift since they 
share the same value Hp = 0, while that of the Roper excitation A^(1440) [up = 1) 
is different. In particular, it is interesting that, although generic values of {mu,md) 
seem to break the isospin invariance, our result shows that the leading order mass 
shift of the baryons is insensitive to this breaking. 

We need to evaluate (p^) with the baryon wave function Rn^ip). The explicit 
expression for the wave function is given in [0. For the lowest Up = 0, it is 



exp 



Mo 2' 



(2.20) 



where Mq = Stt^k, and A/'n^ is the normalization factor. The wave function is nor- 
malized as 

''dpp'R^^{py = l, (2.21) 

with the factor p^ which is the Jacobian for the 4-dimensional spherical coordinate 
system in the one-instanton moduli space. The expectation value is evaluated as 



/OO POO 
dp p'R,{pf / ^ j dt t3/2+2v/TT^g-t 



rdpp'R,{pf V^V rdtfV^^'e-' 
Jo Jo 

_ ( V5Y' V(fe + 3/2) 

"^^'l^2ivj T{b) ' ^ 

where in the second equality we rescaled the variable as p = (-\/6/(2Mo))^/^ t^^^, and 
we have defined 6 = 1 + 2^/1 + N^/5. The numerical value of (|2.22|) for A^^^ = 3 is 
found as 

(P'L=0 /Pel = 2.23. (2.23) 



Using {^^), Mkk = 949 [MeV] and k = aXN^ = 0.00745 which were used in 
2|^, we can evaluate the mass shift ( p.l8| ) for the baryon states with Up = 0. We 



obtain the value of the baryon mass shift as 



1 



ml 



4.11 [GeV^^] . (2.24) 



§This set of values was obtained by fitting /^r = 92.6 [MeV] and the p meson mass in 

^One could use instead the value A/kk — 500 [MeV] which was used in [|l^ for fitting the absolute 

values of the baryon masses, but then it causes an ambiguity on the choice of observable quantities 

one makes for getting the value of A. 
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In particular, our result is universal for all spin, isospin and nz- 

We can compare our result ( |2.24| ) for nucleons {up = 0) with numerical results 
obtained in lattice QCD simulations. For example, a lattice result for nucleons in 
[0 shows ^Mnucieon/"^^ = 4 X 1.02(7) [GeV~^], and a summary of lattice results for 
nucleons in is SM^^cieon/^n = 4 x 0.9 [GeV~^]. Ref. provides a summary 
of lattice QCD results in the year of 2003, and shows values of the coefficient as 
4 X (0.93(5), 1.25(5), 0.93(4), 1.11(4)) [GeV-^] for four different fitting schemes with 
chiral perturbation theory. Rather recent results of lattice QCD [^, ^ show similar 
values. All of these values are consistent with ours ( p.24|) .ll 

A excitation share the same quantum number Hp = 0, thus our value (|2.24| ) also 
applies to A. A lattice result can be found in Ref. as 4 x 1.20(5) for A++ and 
4 X 1.19(8) for A+, which are again consistent with ours ( p. 24 ). 



For Up = 1, the baryon wave function is 



i?^(p) = ATi { _ 1 _ 2^i + iV2/5] p-i+Vi+^eV5exp 



Mo ^ 



(2.25) 



In the same manner, we obtain 



3/2 

3 I V5\ r(6 + 7/2) -26r(6 + 5/2) + 62r(6 + 3/2) 

Pel 



2Nr 



T{b + 2) - 2bT{b + l) + b'^T{b) 



3/2 



3 , 75 \ r(6 + 5/2) + fr(6 + 3/2) 

Pel 



2Nr 



r(6+l) 



(2.26) 



which results in a relation for A^^^ = 3, 



4.16 



(2.27) 



Using this, the mass shift for the baryon states with np = 1 is obtained as 



mt 



7.65 [GeV^^] 



(2.2^ 



II However, note that our result is a large Nc expansion with a sub-leading order term included 
(for distinguishing baryon states), and also obtained under the assumptions which are listed in 



Sec. 2. 1.2. Since the \/Nc corrections are not so small, it is better to think of (2.24) as just an order 
estimate of the baryon mass shift. 
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This is our prediction, stating that the baryon mass shift is larger for the 71^ = 1 
states such as Roper excitation, while for the states with = such as N(1535) 
and A the shift is the same as that of nucleons.** 



3. Discussions on Meson Mass Shift 

The shift of the baryon mass spectrum studied in the previous section relies on an 
assumption that in the large A expansion the term ( p.7|) is the leading order term as 
for the effect of the quark mass term. This term (|2.71) generates a pion mass term as 



seen in Sec. p.l.2| , and the assumption is equivalent to the statement that the shift of 
the vector meson mass spectrum is at a higher order in the large A expansion. In this 
section, we give some arguments for the assumption. The computations given here 
are not sufficient to prove the assumption, but may provide some physical intuition 
on how plausible the assumption is. 

We give two arguments. The first one is based on the chiral perturbation theory 
and its correspondence in string theory amplitudes. The second one is a contribution 
coming from the D6-brane which we introduced to obtain a finite quark mass from 
the worldsheet instanton (see Sec. |2.1.2|) . 



3.1 Chiral perturbation and its stringy realization 



As seen in Sec. p. 1.2 , the string worldsheet instanton induces a term ( p.9| ) of the form 



c tr[Mt/], which is nothing but the term for the quark mass perturbation appearing 
in the chiral perturbation theory. It is the leading term in the derivative expansion. 
In this subsection, we study other possible corrections to the Sakai-Sugimoto effective 
action (p.l|) , from the viewpoint of chiral perturbation theory. 

Sakai-Sugimoto model includes vector mesons as massive gauge fields for hidden 
local symmetries, and then the covariant derivative is derived as 

D^U = d^U - tA^^{x)U + tUA^ix) , (3.1) 

and the gauge fields A^{x) = A{x, z = —oo), A^{x) = A{x, z = +oo) transform as 

Aiix) ^ ^(7^(x)9^(7^(x)-i + (7^(x)A^(7^(x)-i , (?^(x) = g{x, z = -oo) (3.2) 
Aj(x) ^ z(7^(x)5^(?^(x)~i + (7^(x)A^(7^(x)-i , g'^ix) = g{x, z = +oo) , (3.3) 



*It was noted in [g^ that, in the baryon spectrum of massless QCD, there is a degeneracy in 
the mass spectrum for states with common value of the sum Up + nz- In view of our resuhs, this 
degeneracy is resolved once quark masses are turned on. Our results suggest that the mass shift of 
the {np,nz) = (1,0) state (Roper 7V(1440)) is larger than that of the (0,1) state 7V(1535) (which 
is parity odd). However, in reality the observed mass of the Roper is lighter than that of iV(1535). 
This may indicate that the observed resolution of the degeneracy in the mass spectrum is not by 
the quark mass but rather would be by some higher order corrections in 1 /N^ and 1/ A in massless 
holographic QCD. 
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under the gauge transformation g{x, z) G SU{Nf) on the D8-branes. 

Then we can write a higher order term which appears in the sense of chiral 
perturbation theory as 

^TilMUD^U^Df'U] . (3.4) 

This is, in the same manner as tr[M?7]), invariant under the global part of the gauge 
transformation, if we simultaneously rotate the quark mass matrix M, according 
to the philosophy of the chiral perturbation. In the chiral perturbation, this is the 
leading correction to the tr[M?7] term in derivative expansion.''"'' This term in fact 
generates (axial) vector meson mass terms. It is obvious if once we expand the gauge 
field in terms of the KK-decomposed states which are nothing but the vector 
meson states. 

Now, our concern is the order of the coefficient of this term. In the evaluation of 
the baryon mass shift in the previous section, we simply neglected this term. As we 
will see, we can naively argue that this term is at a higher order in the 1/A expansion. 

A stringy interpretation of this term can be found as follows. If we think of this 
( ^.4| ) as a product of MU and D^U^ V^U, the former comes from the worldsheet 
instanton amplitude, while the latter is identical to the pion kinetic term which has 
been obtained in JI], Q from the 5-dimensional YMCS kinetic term. So, we can con- 
clude that this (|3l^) comes from a worldsheet instanton amplitude with fluctuations 
of the gauge flelds on the D8-branes. More precisely, two vertex operators for the 
Yang-Mills fleld on the D8-brane are inserted at the boundary of the worldsheet 
instanton. 

To see the order of this worldsheet instanton amplitude, let us briefly see a 
situation in the flat background. Worldsheet disk amplitudes in flat spacetime are 
of the form 



A — I DX exp I SadXdX-i I da A[X]{dX / da) 

9s J _2tXOl' J /boundary 

The gauge fields here are Taylor-expanded as |^ 



(3.5) 



A[X]{dX/da) = AM[x]{dX^/da) - ^Fmn[x]X^ {dX^ /da) + ■■■ (3.6) 

which is a derivative expansion, and we decomposed X a.s X = x + X where x is the 
zero mode of X[a]. The path integral can be explicitly done to this order because 
the action is quadratic, and the result is nothing but the Born-Infeld action. 



A ^/-det{ri + 27ia'F) . (3.7) 

9s 



ttXhe term Ti[AIDf^D^U] is at the same order, but it is related to (3.4) by the equation 
Udf^Wdf^U ~ d^d^U which can be obtained by a trivial identity 9[C/^C/] = 0. 
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It can be expanded again as 

A^-\l + {a'fF^ + ---] . (3.8) 

9s 

When the disk spans a certain region to give the worldsheet instanton, the zero 
mode X dependence in the semi-classical saddle point in the path-integral generates 
the worldsheet instanton factor exp[— [area]/(27ra')] in front of (p.8|) . 

Now let us try to apply this computation to our case. In our case, we find it 
difficult because the background spacetime is curved, and also because the worldsheet 
boundary is on the D8-branes which are curved and joined. (See (and also PB| ) 
for detailed calculations in a supersymmetric spacetime background.) However, we 
can argue the order of magnitude qualitatively as follows. First, notice that the 
term "1" in (|3.8| ) in fact corresponds to our leading order term ( |2.7|) . The reason 
why we have the path-ordered factor like an Wilson loop in ( p.7| ) is in fact the first 
term on the right hand side of ( |3.6| ). On the other hand, if this ( |3.8| ) is not for a 
worldsheet instanton but a standard disk amplitude {i.e. without winding some D- 
branes), the F"^ term is nothing but D ^U'^ D^U , as has been found in Combining 
these two facts, we come to the conclusion that, compared to the term Mf/, the term 



a 



t\2 



MUD^UW^'U is suppressed by 

Since we are in the warped spacetime, this a' should be replaced by an "effective 
a'" which is given by the tension of a string at at the tip (2; = 0) of the geometry.* 
This is the QCD string tension, 

1 1 



AMj^K- (3.9) 



In the Sakai-Sugimoto model, higher derivative corrections in string theory appears 
with this effective a' (see for example discussions in ||2^). This argument shows that 
the term of our concern is of the form 

cti[MU {l + 0{\-^)D^U^D^'U)] , (3.10) 

corresponding to (|3.8|) . In the large A expansion which the holographic QCD employs, 
the term ii[MU D ^U"^ D^U] is found to be suppressed by A~^. 

This means that (axial) vector meson masses are not so sensitive to the quark 
masses, compared to the pion mass shift. 

3.2 Effects of D6-branes 

In [Q, to regularize the area of the worldsheet instanton, D6-branes ending on the 
flavor D8-branes were introduced. This additional ingredient may affect the meson 
spectra, and in this subsection we discuss the possible effects of the D6-branes. 

*Although the gauge fields on the D8-branes are not only at the tip, the wave functions of 
the vector mesons on the D8-branes is localized at the tip. So we just assume that the effective 
worldsheet propagator in this background geometry may be approximated by the propagator there. 



Note that if one uses a procedure of renormalization of the worldsheet instanton 
proposed in [||, there is no need for the D6-branes, so nothing is worried about con- 
cerning the content of this subsection. On the other hand, the D6-branes in 0] can 
be thought of as a physical cut-off which not only serves as a physical renormaliza- 
tion point but also provides modes on the D6-branes. Therefore, schematically, the 
effective lagrangian of the meson sector in total should be 

'S'dS + 'S'd6 + 5'inst. (3.11) 

In the previous subsection we considered S'inst, and here let us show that the term 
5*06 is sub- leading in 1/A expansion. This is nothing but the D6-brane effective 
action, which is 

Sd6 = -Tm j d'x e- V- ^<9mn + MFmn) ■ (3.12) 

Here the tension of the D6-brane is given by 7d6 = l/(2vr)^/J(7s, and the dilaton in 
the background is e~'^ = (R/U)^^^, with the integral (fx = d^xdTdVL2. Note that 
the D6-brane is along the compactified direction parameterized by r which has a 
period 27r/MKK. 

In the following, we assume that the D6-brane connects the D8 and the anti 
D8-brane straightly without bending. It is expected that the shape of D6 is not 
straight in r but bents, as in the case of the curved shape of the D8-brane. However, 
for simplicity in this paper we assume the straight shape of the D6-brane. 

In addition, in the following computation, we assume that the D6-brane wraps 
the largest S"^ in the S^. As studied in an appendix of [Q, the shape of the D6-brane 
in the U-S'^ space transverse to the N^. D4-branes is non-trivial. However, since our 
purpose here is to estimate the order of magnitude of the D6-brane effect, we may 
take the simple D6-brane shape. 



Expanding the D6-brane action (|3.12|) to quadratic order in the field strength F 



and substituting the background metric and dilaton, we obtain 

^D6 = ^2n)Hlg, j ^'^^^^^2 R"\Uj,,Y'\2^a'f Qf^, + \f{Um)F^^ (3.13) 

where f{U) = 1 — U^^^/U^, and U-qq is the location of the D6-brane in the radial 
direction of the background. Note that the D6-brane is assumed to be straight and 
wrap the 5^, so it is specified by just the coordinate U = Ubg- The action depends 
only on the dimensionless ratio Ubg/Ukk- Introducing a new variable u = U/Ukk, 
we find that the action of the D6-brane, located at u^q = Udq/Ukk, is given by 

^D6 = j^^^ J d'xdrdn, R'l\U^^f'\2r.a'fu]i! Qf^^ + ^(1 - ^5^)^^) . 
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Let us substitute the dictionary for the correspondence between the parameters 
of the supergravity background and the QCD parameters, R = ((^ym^c^s/SMkk)^''^, 
Ukk = {2/9)gYMNcMKKls and = gyuf^'^^KKh- Then we get 

Sue = -^^cMkk I d'xdr u\il (^J^ + ^(l - u^l)F',)l . (3.14) 

Note that all the Is dependence has gone. 

The mode expansion of this action for eigen functions of r gives a spectrum 
of particle modes living on the D6-brane. This is similar to the analysis of flavor 
D-brane fluctuations. On the other hand, since the D6-brane ends on the D8-brane, 
the wave functions of the modes on the D8-branes (which are mesons) can soak into 
the D6-branes, thus have overlap with the D6-brane eigenmodes. We need to look 
at this interaction and diagonalize the quadratic action to see the possible effects on 
the meson mass shift. 

First, let us consider the modes living on the D6-brane. Fourier-decomposing the 
gauge fields along the r direction gives rise to a mass spectrum with a spacing Mkk- 
The eigen functions are trigonometric functions of Mkkt, so the integral / dr in 
( p.l4| ) cancels the factor Mkk in front of the action ( |3.14|) . This means that properly 



normalized fluctuations C^(x) with canonical kinetic terms is 

~ 5^iV-i/2cos(A;MKKr) Cf(a;) (3.15) 

k 

up to a numerical coefficient. The mixing can be evaluated later with this proper 
normalization. 

Next, we look at the meson modes (on the D8-branes) which soak into the D6- 
branes. We have two kinds of mesons: vector mesons whose eigen functions ip2n{z) 
are even in while axial vector mesons whose eigen functions ip2n~i{z) are odd in z. 
Since the D6-branes connect two points on the D8-branes, z = Zm and z = —Zm, so 
depending on the parity of the eigen functions in z, the soaking into the D6-brane 
would be different. For the vector mesons, the most natural soaking would be just 
a constant mode on the D6-brane, while for the axial vector mesons, the natural 
one which costs smallest energy would be cos(Mkkt) on the D6-brane. Once we 
substitute these soaking functions to ( |3.14| ), we obtain 



C / ,4 1/2 

OTT J 



+ liFj-''\z = z^)f + 1m|k(1 - u^l){Af^\z = z^)f 



(3.16) 



Here, we divided the wave functions into the odd part and the even part, 

A^(x,z) = Af'^) + A(r™'^), (3.17) 

oo oo 

^(odd) ^ ^ Bf-\x)ij2n{z), = ^ 5(^-^)(a;)^2n-i(^). (3.18) 



n=l 71=1 
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Note that to obtain the expression ( p.l6|) we identified the gauge fields on the 
D6-brane with that of the D8-branes. The reason why there is no difference in the 
normahzations is as follows. Let us consider a Dp-brane ending on a D(p + 3)-brane. 
The Dp-brane can be thought of as a spike solution of a BPS Dirac monopole on the 
D(p + 3)-brane The (/i = 0, ■ ■ ■ ,p) fiuctuations on the D(p + 3)-brane can 
soak onto the Dp-brane. The DBI action of the D(p + 3)-brane is, if we turn off the 
other fiuctuations, 

Sp+s = rD(p+3) j rf^+"x^det(r/^, + 27Ta'F^,)^det{6ij + diX(-)d,X(-) + 2Tia'Flf) 

where X*^^) and A^^'' are classical solution for the spike, and i, j = j9 + l,p + 2, j9 + 3. 
Explicit substitution of this classical solution gives 

Td(p+3) j dxP^^dxP^^dxP^^^Jdei{5rj + S^X^a^XW + 27ia'Fff) = T^p , (3.19) 



therefore we obtain 



Sp+3 = '^^p j d^^'^x^ det(?7^^ + 27ra' 



(3.20) 



which is the standard Dp-brane action. This shows that the gauge field on the D8- 
brane is in fact a gauge field on the D6-brane, without changing its normalization. 
The normalization of the eigen functions ipniz) is defined in as 



2167r3 

So, defining the normalized eigen function 



dzK ^/''^1pn{z)lpm{z) = Snm- 



we have the standard normalization J dzK~^^^ipn{z)^rniz) 
normalization to the D6-brane action, we obtain 



(3.21) 
(3.22) 

6nm- Substituting this 



^1 



D6 



2^ i'2n-liZm)^2m-liZrn) / d'^X 



A 



n,m=l 



j^p(2n-l)p(2m-l) 



n 2 1/2 



~ /■ 1 1 

^2n{Zm)^2U^m) / d'x -F^'j^^ F^l^^^ + -M^^ - U~^l)B^^-^B^ 

n,m=l 



(2m) 

(3.23) 



Here F are the field strengths of the (axial) vector mesons B^{x). 

The information we need to extract from the last expression (|3.23|) is its overall 
factor 1/A. The action ( p.23|) is supposed to be added to the D8-brane action (which 
is the meson action of [|l|), so we find a 1/A correction due to the D6-brane. 
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There is an overlap between these soaking modes of the mesons and the modes 
on the D6-brane ( p.l5| ). Since the former has a normahzation factor X^^^'^Nc as 
seen m (|]2|) while the latter §J^) has Nc substituting them to we find 



the overlap giving a quadratic interaction of order A"^*^^. 

Therefore, in total, we obtain the kinetic terms for the modes on the D8-branes 
(mesons) and the modes on the D6-brane schematically written as 

Diagonalizing this matrix, we obtain a shift to the vector meson mass which is of 
order 0{\~^^^). As advertised, it is a sub-leading order term, compared to the 
original kinetic terms in the D8-brane action for the vector mesons. Furthermore, 
this shift is further suppressed by the wave function factor ip^Zm)- 



4. Summary 

In this paper, we computed the shift of the baryon mass spectrum by varying the 
quark mass, from the chiral limit, in the Sakai-Sugimoto model 0] of holographic 
QCD. To introduce the quark mass to the model, we employed the worldsheet in- 
stanton approach 0, ^. We combined it with the solitonic description of the baryon 
used in [jlO[, and we obtained analytically the mass shift of the baryons, which is 



linear in the pion mass squared. The slope depends on the species of the baryons. 

The slope computed for nucleons (whose value is also shared with A), with inputs 
chosen to be and nip as in Jl], 0, roughly agrees with lattice data for nucleons 
and A. We find that the slope is shared also with A^(1535), while the slope for the 
Roper excitation A^(1440) is found to be larger. These slopes for excited baryons 
are our predictions, and it would be quite interesting if they can be confirmed by 
lattice computations in the future. Although we work in large Nc QCD with the 
holographic approach, we hope that the tendency of the difference in the slope may 
be the same for QCD. 

In Sec. 3, we gave various discussions on the vector meson mass shift, but our 
naive arguments showed that possible effects are subleading order in the large A 
expansion (which justifies our computation in Sec. 2). It would be important to 
clarify this issue further and compute the vector meson mass shift more rigorously. 
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